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It is well known that there are resonant states with complex energy for the supercritical Coulomb 
impurity in graphene. We show that opening of a quasiparticle gap decreases the imaginary part 
of energy, |Im_E|, of these states and stabilizes the system. For gapless quasiparticles with strong 
Coulomb interaction in graphene, we solve the Bethe-Salpeter equation for the electron - hole bound 
state and show that it has a tachyonic solution for strong enough coupling a = e 2 /kHvf leading to 
instability of the system. In the random phase approximation, the critical coupling is estimated to 
be a c = 1.62 and is an analogue of the critical charge in the Coulomb center problem. We argue 
that the excitonic instability should be resolved through the formation of an excitonic condensate 
and gap generation in the quasiparticle spectrum. 



I. INTRODUCTION 

Graphene is an one-atom-thick layer of graphite packed in the honeycomb lattice. Although theoretically considered 
long time ago [l[ , graphene became an active area of research only recently after the experimental fabrication Q of 
this material and because of a variety of its unusual electronic properties. 

At low energy the band structure of graphene is formed by the 7r-electron orbits of carbon and consists of the valence 
and conduction bands and corresponds to cones touching each other at the so-called Dirac points. Quasiparticle 
excitations close to these points are described by the massless Dirac equation and have a relativistic-like dispersion 
E = ±fi,«i?|k|, where vf ~ 10 6 m/s is the Fermi velocity and k is the quasiparticle wave vector. This fact brings an 
exciting connection between graphene and 3 + 1- dimensional quantum electrodynamics (QED). 

The vanishing density of states at the Dirac points ensures that the Coulomb interaction between the electrons in 
graphene retains its long-range character in view of vanishing of the static polarization function for q — > Q. The 
large value of the coupling constant a = e 2 jfivp ~ 1 means that a strong attraction takes place between electrons and 
holes in graphene and this resembles strongly coupled QED, thus providing an opportunity for studying the strong 
coupling phase experimentally within a condensed matter laboratory. Given the strong attraction, one may expect 
an instability in the excitonic channel in graphene with subsequent quantum phase transition to a phase with gapped 
quasiparticles that may turn graphene into an insulator. This semimetal-insulator transition in graphene is widely 
discussed now in the literature [JQ since the first study of the problem in Refs.0, 0|- The gap opening is similar to 
the chiral symmetry breaking phenomenon that occurs in strongly coupled QED and was studied in the 70-tics and 
80-ties @, H, E3, EH; E3 ■ I n fact, the predicted strong coupling phase of QED, like other QED effects not yet observed 
in nature (Klein tunneling, Schwingcr effect, etc.), has a chance to be tested in graphene. 

We begin our study with the problem of the supercritical Coulomb center in Sec. [II] m graphene. As is well 
known [la, fl4|. for the Coulomb potential, Vc(r) = —Ze 2 /nr, the spectrum of quasiparticles with a gap A contains 
a continuum spectrum for \E\ > A and a discrete one for < E < A. The lowest bound state energy equals 
Eq = A-^/l — (2Za) 2 and becomes purely imaginary for Za > 1/2 - the "fall into the center" phenomenon. The 
unphysical complex energies indicate that the Hamiltonian of the system is not a self-adjoint operator for supercritical 
values Za > 1/2 and should be extended to become a self-adjoint operator. The way out of this situation is well 
known from the study of the Dirac equation in QED: one should replace the singular 1/r potential by a regularized 
potential which takes into account the finite size of the nucleus, i?, [HI, [HI, [13] ■ When the charge Z increases, the 
energies of discrete states approach the negative energy continuum, E = —A, and then dive into it. Then discrete 
states turn into resonances with a finite lifetime, which can be described as quasistationary states with complex 
energies, ImE ^ 0. Such states correspond to a rearrangement process when an electron-hole pair is created from the 
vacuum, the positively charged hole goes to infinity and the electron is coupled to the Coulomb center, thus shielding 
the charge of the latter. The critical charge Z c is determined by the condition of appearance of nonzero imaginary 
part of the energy and increases with the increase of A. 

Turning to the case of gapless quasiparticles in the regularized Coulomb potential, there are no discrete levels for 
Za < 1/2 due to scale invariance of the massless Dirac equation, and for Za > 1/2 quasistationary states emerge [UJ] ■ 
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The energy of quasistationary levels for the regularized potential has a characteristic exponential-type dependence, 
HeE, ImE ~ — R^ 1 cxp(— ir/^/(Za) 2 — 1/4), in the nearcritical regime (according to the analysis in Appendix A, the 
critical coupling Z c a — > 1/2 for RA — > 0). We find that switching on a fcrmion gap, A <C \E\, decreases |Imi?|, i.e. 
increases the stability of the system. The situation here is analogous to the problem of a massless electron in the 
supercritical Coulomb center in QED first studied in (for a review, see, @). 

In Sec lIIII we show that the instability in the supercritical Coulomb center problem is closely related to the excitonic 
instability in graphene in the supercritical coupling constant regime a > a c ~ 1. Solving the Bethe-Salpeter equation 
for an electron - hole bound state in graphene, we demonstrate that for strong enough coupling constant there are 
tachyon states with imaginary energy (E 2 < 0) in the spectrum which play here the role of the quasistationary states 
in the problem of the supercritical Coulomb center. The presence of tachyons signals that the normal state of freely 
standing graphene is unstable. In fact, the tachyon instability can be viewed as the field theory analogue of the "fall 
into the center" phenomenon and the critical coupling a c is an analogue of the critical coupling constant Z c e 2 /hvp in 
the problem of the Coulomb center. However, in view of the many-body character of the problem, the way of curing 
the instability in graphene (like in QED [§]) is quite different from that in the case of the supercritical Coulomb center. 
Since the coupling constant in freely standing graphene a ks 2.19 is larger than 1/2, the quasielectron in graphene has 
the supercritical Coulomb charge. This leads to the production of an electron-hole pair, the hole is coupled to the 
initial quasielectron forming a bound state but the emitted quasielectron has again a supercritical charge. Thus the 
process of creating pairs continues leading to the formation of exciton (chiral) condensate in the stable phase, and, as 
a result, the quasiparticles acquire a gap. The exciton condensate formation resolves the problem of instability, hence 
a gap generation should take place in a freely standing graphene making it an insulator. 

Sec. IV contains our conclusions. In Appendix [5] we consider the behavior of bound states for gapped graphene 
quasiparticles in the regularized Coulomb center and find the critical coupling Z c a as a function of the parameter RA 
for the lowest energy level. In Appendix [B] we give the exact solution for the tachyon wave function which satisfies 
the fourth order differential equation. 



II. THE SUPERCRITICAL COULOMB CENTER: RESONANT STATES 



Although the electron-hole problem is a many body problem in graphene, it is instructive to consider a rather 
simple one particle problem of the electron in the field of the supercritical Coulomb center in view of the connection 
of the latter problem with the excitonic instability in graphene (for a similar problem of instability in the case of a 
massless fermion in the external field of the supercritical Coulomb center in QED see Refs. [9jjl9j) . 

The Coulomb center problem was studied quite in detail in the literature [HI, [U, [H, l20l. l2ll. 1 23. HH, [24| ■ Since we 
are mainly interested here in resonant states, we will consider a regularized Coulomb potential 

V(r) = - — , (r > R), V(r) = (r < R), (2.1) 

nr kR 

(k is a dielectric constant) because resonant states are connected with diving into the lower continuum that takes 
place only in the case of a regularized potential [l7], [H[ . 

The electron quasiparticle states in vicinity of the K point of graphene in the field of Coulomb impurity arc described 
by the Dirac Hamiltonian in 2+1 dimensions 

H = (cr 3 A + V(r) - ihv F a x d x - ihv F a 2 d y ) , (2.2) 

where a 1 are Pauli matrices, vp the Fermi velocity. (The Hamiltonian of quasiparticle excitations near the K' point is 
given by (|2.2[) with matrices a 1 multiplied by -1.) Note that we introduced the Dirac gap A. Although it is absent in 
the quasiparticle Hamiltonian in graphene in view of the E/(4) symmetry, it may appear due to spontaneous symmetry 
breaking. Since Hamiltonian (|2.2[) commutes with the total angular momentum operator J z = L z + S z = — *^^ + fc3, 
we seek eigenfunctions in the following form: 

1 / e #(i-V2) a ( r ) 



^» e ^0-+i/3)6( r ) ) ■ ( 2 - 3 ) 
Then we obtain a system of two coupled ordinary differential equations of the first order 

s-u+imi+ z+tzmt -o, v+ij - m i- ^=m a=0 . (2.4) 

r Hvp r hvp 
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It is convenient to define the variables e = E/Hvf, m = A/Hvf, 



\prrP- 



2ur, and a — e 2 jtmpK. Equations 



are solved in Appendix [Al where the discrete spectrum is found in the weak coupling regime Z < Z c . According 
to the analysis there, the critical coupling Z c a — > 1/2 for mR — ► 0. 

Let us analyze Eqs. ([2.4[ ) in the supercritical case Za > 1/2 and show that there are resonant states for |e| > m 
(we define the gap A > 0). These states describe the instability of the supercritical charge problem with respect to 
the creation of electron-hole pairs from the vacuum. The created electron is coupled to the Coulomb center, thus 
shielding the charge of the latter while positively charged hole goes to infinity [l?], HH ; the process is repeated until 
the charge of the Coulomb center is reduced to a subcritical value. 

The Whittakcr function W fJ ,.^(p) with p = 1/2 + Zae/u, v = y^j 2 — Z 2 a 2 describes bound states for |e| < |m| which 



are situated on the first physical sheet of the variable u and for which Reu > (see, Eq. (|A7[0 . The quasistationary 
states are described by the same function W fliU (p) and are on the second unphysical sheet with Rett < 0. We shall 
look for the solutions corresponding to the quasistationary states which define outgoing hole waves at r — > oo with 



Ree < 0, line < 0, Rcu < 0, Imu < 0. 



(2.5) 



For solutions with Z 2 a 2 > j 2 resonance states are determined by Eq. (|All[ ) for bound states where v is replaced by 
v = i/3. We will consider the states with j = 1/2 which correspond to the nSi -states, in particular, the lowest energy 
state belongs to them. The corresponding equation then takes the form 



A3 



(p) 



\2 u j 



w 



r=R 



fc+1 



Za/ R — m Ji(p) 
e + Za/ R + m Jo(p) 



y/(Za + eR) 2 - m 2 R 2 , (2.6) 



where W^ v (x) and J a (x) are the Whittaker and Bcsscl functions, respectively 

We are interested in the case of |e| 3> m, and more important, in the case of the massless electron, m = 0. The 
analytical results can be obtained for the near-critical values of Z when Za — 1/2 -C 1. We assume that \2uR\ <C 1, 
then using the asymptotic of the Whittaker function, we find 



(2uR) 



2hi 



r(l - 2i/3) r (1 + i/3 



i/3 



Za(r 



e) I 
— 2 



i/3 — Za 



Jl(Za) 
Jo(Za) 



T(l + 20) T{1- i/3 I 
Expanding Eq. (|2.7[) in the near critical region in powers of (3 



i/3 



Za(r 



!) 1 

2 



i(3 



Jo(Za) 



(2.7) 



yf Z 2 a 2 — 1/4, we find the following equation: 



(-2iy/e 2 - m 2 R) 210 = 1 + 40 



Jo(l/2) 



J (l/2)- Jr(l/2) 



*(1) 



1 



2 Ve 



A/ e+m 



(2. 



Here *f?(x) is the psi-function and we put u = e 2 — m 2 where Im-\/ e 2 — m 2 < on the second sheet. At first we 
consider the case m = 0. Writing e = \e\e n Eq. (|2.8|) takes the form 



ln(2\e\R) + i[ 7 --)~2 



Jo(l/2) 



Ml/2)- Ml/2) 



* ( i) _ f i _ ^ _ JL 

V ' 2 \ 2 1 + i 



nn 

- T , n=l,2, 



(2.9) 



We find 



,(0) 



airV 7 exp 



v/Z 2 a 2 -l/4 



(1.18 + 0.17z)i? _1 exp 



y/Z 2 a 2 - 1/4 



n= 1,2,.. 



(2.10) 



where 



7T / 7T 

7 = — 1 + coth — 
' 2 V 2 



1 

a = — exp 
2 



3.28, 
2J (l/2) 



N + 2*(1) - 1 - Re* 1 
Jo(l/2) - Ml/2) 1 j V 



1.19. 



(2.11) 
(2.12) 



The energy of quasistationary states (|2.10[) has a characteristic essential-singularity type dependence on the coupling 
constant reflecting the scale invariance of the Coulomb potential. The infinite number of quasistationary levels is 
related to the long-range character of the Coulomb potential. Note that a similar dependence takes place in the 
supercritical Coulomb center problem in QED [lj|. Our results are also in agreement with Ref . fl8l] . 



4 



Since the "fine structure constant" e 2 /Hvf ~ 2.19 in graphene, an instability appears already for the charge Z = 1. 
However, in the analysis above we did not take into account the vacuum polarization effects. Considering these effects 
and treating the electron-electron interaction in the Hartree approximation, it was shown in Ref. that the effective 
charge of impurity Z e ff is such that the impurity with bare charge Z = 1 remains subcritical, Z e f fe 2 / (nhvp) < 1/2, 
for any coupling e 2 /(kHvf), while impurities with higher Z may become supercritical. 

For finite m and in the case |e| 3> m, Ree < 0, expanding Eq. (|2.8|) in m/e we get up to the terms of order m 2 /e 2 , 



m 
~2e 



(?) i _ 



-(0.29 - 0.23i) 



n = 1,2, 



(2.13) 



The resonant states with describe the spontaneous emission of positively charged holes when electron bound states 
dive into the lower continuum in the case m = 0. In order to find corrections to these energy levels due to nonzero 

oo 

k and easily find the first two terms 



m, we seek solution of Eq. (|2.13p as a series e = with of order rr 

k=0 



& - m 



(0.24 + 0.2CK). 



(2.14) 



Since Ime { n ] < 0, the appearance of a gap results in decreasing the width of resonance and, therefore, increases 
stability of the system. 

We considered above the case |e| 3> m and analyzed how a nonzero mass affects resonant states. It is instructive to 
consider resonant states also in the vicinity of the level e = — m when bound states dive into the lower continuum and 
determine their real and imaginary parts of energy. First of all, nonzero m increases the value of the critical charge. 
Let us find it. Using Eq. (|A20|) in Appendix A, we obtain that the critical value Z c a for j = 1/2 scales with m like 
(see Fig.l) 



Z r a 



1 



2 log 2 (cmR) 



c = exp 



-2^(1) - 



2J (l/2) 



Jo(l/2)-Ji(l/2) 



0.21. 



(2.15) 



Note that the dependence of the critical coupling on mR is quite similar to that in the strongly coupled QED [1, 0] . 
For Z > Z c , using Eq. (|2.8|) . we find the following resonant states: 




FIG. 1: The critical coupling as a function of mR for the IS1/2 level. 



where f3 c = *J (Z c a) 2 — 1/4. Like in QED [l6[ the imaginary part of energy of these resonant states vanishes ex- 
ponentially as Z — > Z c . Such a behavior is connected with tunneling through the Coulomb barrier in the problem 
under consideration. For the quasielectron in graphene in a central potential V(r), expressing the lower component 
of the Dirac spinor (|2.3p through the upper one and following [l6|, [l7| , we obtain an effective second order differential 
equation in the form of the Schrodinger equation 

x"( r ) + k 2 (r)x(r) = 0, a(r) = exp 



V 



2J y r e + m 



V 



I dr 



X{r) 



(2.17) 
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Here 



k 2 {r) =2(£ -U(r)), £ 



V 

Hvf 



(2.18) 



and we represent the effective potential as the sum of two terms U = U\ + U2, where U\ is the effective potential for 
the Klein-Gordon equation and Ui takes into account the spin dependent effects 



- V 2 
U X = eV-- 



2r 2 : 



U 2 



V" 



V 



2jV 



m — V 2 \ e + m — V 



-(e + m - V) 



(2.19) 
(2.20) 



Note that Eq. (|2.17p and the potentials (|2. 19[) , (|2.20|) coincide with the corresponding equations in QED [l7| • We plot 
the effective potential U(r) for Z — > Z c , j — 1/2, and e = —m in FigfJl where the Coulomb barrier is clearly seen. 




FIG. 2: Effective potential for the Coulomb center in the case e — —m and Z — Z c 



Up to now, we considered a one particle problem in an external field. In the next section, we will consider electrons 
and holes in graphene which interact by means of the Coulomb interaction and show that an instability develops in 
the system when the coupling a exceeds some critical value a c . 



III. THE EXCITONIC TACHYON INSTABILITY 



A. The Bethe— Salpeter equation 



The instability of the supercritical charge problem due to the emission of positively charged holes discussed in the 
previous section indicates the possibility of the excitonic instability in graphene in the case of a supercritical coupling 
constant. In this section, we will study the Bethe-Salpeter (BS) equation for an electron-hole state and show that it has 
a tachyon in the spectrum in the supercritical regime. Before wc do this, let us discuss some similarities and differences 
of the supercritical Coulomb charge problem with the famous Cooper problem in the theory of superconductivity. 

Although the Cooper problem is formulated as a quantum mechanical problem for two particles (electrons), it can 
be standardly reduced to a one particle problem in an external potential. Therefore, the Coulomb center problem is 
similar to the Cooper problem in this respect. However, there are important differences between the two problems. 
The first one is connected with the fact that the Dirac equation contains the lower continuum with filled negative 
energy states. Therefore, if a bound state energy enters the lower continuum, we are essentially dealing with a many 
body problem. This explains why there are resonant states with imaginary energy in the supercritical Coulomb 
potential unlike the Cooper problem where there are only negative energy bound states. The second important 
difference between these two problems is connected with the critical value of coupling constant. It is zero for the 
Cooper problem because the density of states in this problem is nonzero at the Fermi surface that plays a crucial role 
in the bound states formation. On the other hand, a c = 1/2 for the Coulomb center problem in graphene where the 
density of states is zero at the Dirac point. 

The appearance of the Cooper bound state in the theory of superconductivity is directly related to the instability 
of the normal state of metal. Indeed, according to [2(|, the BS equation for an electron-electron bound state in the 
normal state of metal has a solution with imaginary energy, i.e. a tachyon. This means that normal state is unstable 
and a phase transition to the superconducting state takes place. As we mentioned above, resonant states in the 
supercritical Coulomb center problem suggest the excitonic instability in graphene. 
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For the description of the dynamics in graphene, we will use the same model as in Refs.@, 0] in which while 
quasiparticles are confined to a two-dimensional plane, the electromagnetic (Coulomb) interaction between them is 
three-dimensional in nature. The low-energy quasiparticles excitations in graphene arc conveniently described in 
terms of a four-component Dirac spinor ty^ = {tl>KAai'4 , KBai'4 ) K'Bai' i pK'Aa) which combines the Bloch states with 
spin indices a = 1, 2 on the two different sublattices (A, B) of the hexagonal graphene lattice and with momenta near 
the two nonequivalent valley points (K, K') of the two-dimensional Brillouin zone. In what follows we treat the spin 
index as a "flavor" index with Nf components, a = 1, 2, . . . Nf, then Nf = 2 corresponds to graphene monolayer while 
Nf = 4 is related to the case of two decoupled graphene layers, interacting solely via the Coulomb interaction. 

The action describing graphene quasiparticles interacting through the Coulomb potential has the form 



S = J dtd 2 r^ a (t, r) (i^dt - iv Fl V) * a (t, r) 

- i y"d^ C ^<'d 2 rd 2 /^ a (^,r) 7 * a (^,r)[/ (^-^^|r-r'|)* b (^',r')7 * 6 (^^r'), (3.1) 

where VP = and the 4x4 Dirac 7-matrices 7^ = t 3 (a 3 , icr 2 , — ier 1 ) furnish a reducible representation of the 

Dirac algebra in 2 + 1 dimensions. The Pauli matrices r, a act in the subspaces of the valleys (K, K') and sublattices 
(A, B), respectively. The other two 7-matrices which we use are 7 3 = ir 2 <g> (Jo, 7 5 = i7°7 1 7 2 7 3 = t\ <E> cr (<r is the 
2x2 unit matrix). 

The bare Coulomb potential Uo(t, |r|) takes the simple form: 

However, the polarization effects considerably modify this bare Coulomb potential and the interaction will be 

IJU Irh - e ' [ dLJ ( rf2fc cxp(-^ + ikr) 

^^ |r|j " K j2nJ 2tt |k|+n(w,k) ' ^ 

where k is the dielectric constant due to a substrate and the polarization function II (a;, k) is proportional (within the 
factor 2tt/k) to the time component of the photon polarization function. Correspondingly, the Coulomb propagator 
has the form 



The one-loop polarization function is Q 



D ^=wrm^wr (3 - 4) 



. , , ne 2 N f k 2 . , 

4k y/h 2 v 2 pk 2 - u} 2 



and in an instantaneous approximation it is 



In general, the static polarization operator must have the form 11(0, |q|) = \q\F(a, Nf) due to dimensional reasons, 
however its exact form is not known and in the present paper we will use the one-loop approximation. 

The continuum effective theory described by the action (|3.ip possesses the U(2Nf) symmetry. However, as was 
pointed out in Ref. 27 1 (see also Refs. [H, H^j), it is not exact for the Lagrangian on the graphene lattice. In fact, 



there are small on-site repulsion interaction terms which break the U(2Nf) symmetry. 

In order to analyze excitonic instability, we consider the Bethc-Salpcter equation (see, for example, Ref. [12||) for an 
electron-hole bound state which is represented in Fig[3] The kernel K of the BS equation in the simplest approximation 
contains two diagrams: the one is due to exchange Coulomb forces and another one is the annihilation diagram. The 
annihilation diagram does not contribute for the BS wave function considered below. Thus the BS equation takes the 
following form: 



S- 1 {q+ l -P) X (q,P)S- 1 {q- l -P) 



a/3 



yd 3 fcD(|q-k|)[7 x(fc,P)7°] a/3J (3.7) 





a: 



— J — 1 




1 — > — 


— <— 




— <— 




FIG. 3: The BS equation for a bound electron-hole state \. The kernel if contains two diagrams: exchange and annihilation 
ones. The wave line corresponds to the Coulomb propagator. 

where k = (fco,k), a,j3 are spinor indices, xilyP) is the BS amplitude in momentum space 



\P), 



(3.8) 



q = (qo,ci), P = (Po,P), q and P are relative and total momenta, respectively, and 

■A 



S(P) 



TPo ~ 7P - 



p 2 -p 2 -A 2 + iO 



is the quasiparticle propagator with a gap A (the gap A is zero in non-interacting graphene, however, may be generated 
due to the strong Coulomb interaction). In what follows we put h = vp = 1. 

Taking into account the static vacuum polarization by massless fermions, i.e., with LT(c<j = 0,k), corresponds to the 
replacement of the coupling constant a in Eq. (|3.7p by 



Further, introducing the function 



1 + 7ra7V//4 



1 



2A. 



x{q,p) = S-\q+-P)x{q,P)S- 1 {q 



-P) 



with p = P/2, the BS equation can be equivalently rewritten as follows 

2iX f d 3 k 



x(q,p) 



1 °S(k + p)x(k,P)S(k-p)y°. 



(3.9) 



In general, x can be expanded in 16 independent matrix structures. In view of the experience in QED 0], we expect 
a gap generation in graphene in the supercritical regime. Then the spin- valley U (4) symmetry will be broken (see, 
eg- [a 0]) that leads to the appearance of massless Nambu-Goldstone bosons in the spectrum. Similarly to QED 
[9j , these Nambu-Goldstone bosons are transformed into tachyons if considered on the wrong vacuum state without 
a gap generation. In the present paper, we will consider only matrix structures of x connected with the 7 5 matrix 



x(q,p) = X5{q,P)l b + X05(q , ,p)9 i 7 i 7°7 5 



(3.10) 



where Xb{qiP) an d Xo5(q,p) are scalar coefficient functions. We will see in the next section that it is enough to consider 
only Xb m order to describe a Nambu-Goldstone excitation in the massive state. However, we retain the function X05 
because it is necessary in the study of tachyon. There are also tachyons in other channels which describe different 
ways of breaking the U(2Nf) symmetry, for example, one can use matrices /,7 3 ,7 3 7 5 instead of the matrix 7 5 in 
Eq. (|3.10[) . To study instability it is enough to find at least one channel with tachyons. The real pattern of a symmetry 
breaking is defined by solving gap equations for various kinds of order parameters and determining which of them 
corresponds to the global minimum of the system energy. For simplicity we consider only the channel described by 
the wave function (|3.10|) which can be treated analytically. 
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B. Tachyon states 



Let us first show that, for A > A c , there is a tachyon in the spectrum of the Bethe-Salpeter equation in the 
masslcss theory A = and determine the critical value A c . For the study of tachyon, we can set p = 0, however, 
should keep nonzero po- One can check that ansatz (|3. 10|) is consistent for Eq. (|3.9p and leads to a coupled system of 
equations for functions Xb{qiPo)i Xo5(q,Po) (in what follows we omit po for brevity in the arguments of the functions 
Xb{q->Po), XosfePo))- Since Eq. (|3.9p implies that x(q,p) does not depend on q , we can integrate then over k by using 
the integrals 



dko 



ci +c 2 k + c 3 k$ 



* ((*b - Po) 2 - k 2 + iS) ((fc + Po) 2 - k 2 + id) 

— oo 

where 5 — > +0. We obtain the following system of integral equations: 



ci + c 3 (p§-k 2 ) 
2|k|( P 2 -k 2 ) ' 



xs(q) = A 



d 2 fck 2 ( X5 (k)+p X05(k)) 



xos(q) = A 



d 2 fcqk(k 2 xo 5 (k) +p x 5 (k)) 
2tt q 2 |q-k||k|(k 2 -p 2 ) 



27T |q-k||k|(k 2 -p 2 ) ' 
We assume that Xs(q) and xos(q) depend only on q = |q|, then we can integrate over the angle using 

2tt 



(1,0 





/ 



1/g 2 + fc 2 — 2qk cos < 
cos </> 

■\/g 2 + fc 2 — 2qfc cos < 



-K 



2\/qk 
q + k 



9(q - k) 



K 



<1 



e{k-q) K (q 



(3.11) 



(3.12) 



2(g 2 + fc 2 ) 
qk(q + k) 

q0(q - k) 



2v^fc\ (q + k) 2 ^/2^k 



4 

qk 



q + k 
k 



k 2 



E 



- E 



(3.13) 



where if (a:) and are complete elliptic integrals of the first and second kind, respectively, 9{x) is the Heaviside 

step function, and for the last equalities in Eqs. (|3.12p . (|3. 13(1 we used the formulae 8.126.3,8.126.4 in the book [30j . 
Further, we approximate the elliptic integrals by their asymptotics 



%)^(l + 0(x 2 )), E(x)~^(l + 0(x 2 )), x«l. 



(3.14) 



This approximation allows one to obtain analytical results for the BS equation. The logarithmic singularity present 
in elliptic integrals in Eqs. (|3.12p . (|3. 13|) at q = k does not influence qualitatively the solution obtained though it is 
important to take it into account to get correct value of the critical coupling (see, the derivation of Eq. (|3.30p below). 
Thus we find 



X5(q) = A 



k 2 dk 



q(k 2 -p 2 ) 



(Xs(k) +PoXo5(k)) + A 



kdk 

■ 2 -pI 



(X5(k) +PoXo5(k)) , 



(3.15) 



xo5{q) 



k 2 dk 



2 J q*{k 2 -p 2 ) 




^ (k 2 X05(k) +PoX5{k)) + - 



dk 



k(k 2 -p 2 ) 



{k 2 Xoz{k) +poX5(k)) 



(3.16) 



Here we also introduced a finite ultraviolet cutoff A which could be taken to be of order tt / a, where a is a characteristic 
lattice size, a = 2.46A for graphcnc. An alternative, equally good, choice of A is related to the energy band, A = t/vp, 
where t = 2.4 eV in graphene. 

These equations are equivalent to the system of differential equations 



// . 2 , X5 + P0X05 

X 5 + -Xb + A — 2 — 

q q - Po 



o. 



with the following boundary conditions: 

0, [qX5(q))' 



q\' 5 



q=0 



,, , 4 , 3A g 2 XQ5 + PoX5 

X05 ~ ' X05 ' T, 2~T~2 TT 

q 2 q 2 {q z -pt ) ) 

= 0, [q 3 Xo 5 (q)Y = 

q=0 q=A 







0, <7 4 Xo5 



(3.17) 



(3.18) 
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The system of differential equations (|3.17|) can be reduced to one equation of the fourth order whose solutions are 
given in terms of generalized hypergeometric functions iF^(q 2 /p 2 ,) and the Meijer functions with the corresponding 
boundary conditions (for this analysis see appendix |B|) . However, since we seek for the solution with po — > 0, it is 
simpler to analyze straightforwardly the system (|3.17p itself, in this regime the system decouples 



X5 + 7X5 + A- 



A'5 



Po 



0, 



X05 + -X05 + 



3A_Xo5 
2 q 2 



0. 



(3.19) 



where we keep po in the denominators because it regularizes singularities for q — > 0. 

Obviously, Eqs. (|3.19p are differential equations for the hypergeometric function F(a, b; c; z) [30| . The solutions that 
satisfy the infrared boundary conditions are 



T + 7 1-7 3 q 2 \ ^ 3(1 + 7) 3(1-7) 5 q 2 

XB = CiF[- r ,—,- l ^J, Xo5 = C 2 F^-^,^^;-,^ 



(3.20) 



where 7 = yl — 4A and 7 = y/l — 2A/3. Using the asymptotic of the hypergeometric functions, one may easily check 
that the ultraviolet boundary conditions for the function \5 can be satisfied only for A > 1/4, therefore, 1/4 is the 
critical coupling for the approximation that we use. (Note that if we neglect the vacuum polarization contribution, 
then A = a/2 and the critical value 1/4 coincides with the critical coupling Z c a = 1/2 obtained in Sec. II for the 
Coulomb center problem.) The UV boundary condition for the function X05 can be satisfied for the values of A > 3/2 
but not for A < 3/2. Therefore, for 1/4 < A < 3/2 we take a trivial solution X05 = and we are left only with the 
equation for the function \5- Knowing the function \5 wc then solve an inhomogencous equation (|3.17jl for X05, in 
this way we find that the function X05 ~ Po ■ The critical value A c = 1/4 coincides with the critical coupling constant 
found in Q, where the same approximation for the kernel was made. In the supercritical regime 7 = iu), to = y/4X — 1 
and the function xs(<7) behaves asymptotically as 



X5{q) 



q 1 / 2 cos 



(\A - 1/4 In q + constj 



(3.21) 



Such oscillatory behavior is typical for the phenomenon known in quantum mechanics as the collapse ("fall into the 
center") phenomenon: in this case the energy of a system is unbounded from below and there is no ground state. 
Nodes of the wave function of the bound state signify the existence of the tachyon states with imaginary energy po, 
Imj»Q < 0. Indeed, the UV boundary condition for \5 leads to the equation 



(i + iw)r(i + f )r 



' l-iu 



r 



' 5 — iuj N 



(1 - iu)T (1 - f ) r (±±p) r (£±s£) 

Then wc find the following tachyon solution: 



4 1 

Po 



(3.22) 



47m 

Po = —A cxp h S(u) 



8(u>) = - 



If A tends to 1/4 from the above, i.e. ui — > 0, 



arctanw + Arg ( F ( 1 + — ) T 



l-iu 



5 — iu) 



(3.23) 



Pa - = -A z atp[ — — + <5(0)), 5(0) =4 + 2*(l)-*(l/4)-*(5/4) w 7.3, n = l,2,.... (3.24) 



Thus, we see that the strongest instability, i.e., the smallest negative value of p\ is given by the solution for the function 
Xb with n = 1. The tachyon states play here the role of the quasistationary states in the problem of supercritical 
Coulomb center resulting in the vacuum instability. In fact, the tachyon instability can be viewed as the field theory 
analogue of the "fall into the center" phenomenon and the critical coupling a c is an analogue of the critical coupling 
Z c a in the problem of the Coulomb center. 

The tachyon energy Eq. (|3.24[) has a characteristic essential singularity of the kind 1/\A — A c in the exponent. It 
can be argued that this behavior reflects a scale invariance in the problem under consideration and keeps its form for 
any approximation which does not introduce new scale parameter except the cutoff [31 1. 

There are two possibilities for the system with the supercritical charge to become stable: to shield spontaneously 
the charge or to generate spontaneously the fermion gap. The first possibility is realized in the problem of the 
supercritical Coulomb center which is due to the formulation of the problem as the one-particle one. The second 
possibility - dynamical generation of the fermion gap - is realized for quasiparticles in graphene interacting through 
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supercritical Coulomb interaction. The situation here is completely analogous to the strongly coupled QED [1, 0, [l2j 
where it is shown that the vacuum stabilization by generating dynamical fermion gap is a rather universal phenomenon. 
The critical value A c determines the critical coupling a c as a function of the fermion number Nf, 



4A n 



(3.25) 



2-TriV/Ac' 

(compare with Eq.(28) in Q). The critical value A c = 1/4 in the approximation (|3. 14)) used for kernels. The more 
precise value of A c can be found if one notes that A c corresponds to the limit po = 0. Taking this limit in the system 
(j3~T7|) . we get 



X5{q) 



2A 

— / dkxn{k) 



TT 



9(q - k) 



K I - 

q 



9{k - q) 



K 



(!) 



(3.26) 



Xon(q) 



2A 
71 ~q 



dkxob{k) 



9(q -k)[K 



- E 



ke(k-q) 



- E 



(3.27) 



Note that the ultraviolet cutoff, A, has been taken to infinity, which is appropriate at the critical point. These 



equations are scale invariant and are solved by Xs(q) 
satisfy the transcendental equations 



q J iXo5(q) = q p on the condition that the exponents j,p 



— f dx\a 

TT J 

o 

l 



+ x~<- 1 \K(x), 0< 7 <1, 



x p - 3 \ (K{x) - E(xj), < p < 3. 



(3.28) 



(3.29) 



These equations define roots 7, p for any value of the coupling A. An instability is signalized by oscillatory behavior of 
the functions X5(<7)> Xosiq)- F° r the function Xb{q) this occurs when two of the roots of Eq. (|3.28p in the interval (0, 1) 
coalesce and then become complex conjugate. We find that this happens when 7 = 1/2, for this value the integral in 
Eq. (|3.28p is exactly evaluated (see, the book [32|) and we obtain the critical value 



A c = 



4tt 2 



r 4 (i/4) 



0.23. 



(3.30) 



The second equation (|3.29p gives higher critical value A c = 0.91 therefore the instability is determined by the value 
A c = 0.23. The critical value N cr i t 2.8 corresponds to a = 00 in Eq. (|3.25p . Since for graphene the number of 
"flavors" Nf = 2, the critical coupling is estimated to be a c w 1.62 in the considered approximation (3?]]. Because 
the coupling constant in freely standing graphene a ~ 2.19 (k ~ 1) the system is in the unstable phase. On the other 
hand, for graphene on a SiC>2 substrate the dielectric constant k ss 2.8, therefore, a f» 0.78, i.e., the system is in the 
stable phase. 

Finally, since the U(2Nf) symmetry is spontaneously broken, there must exist Nambu-Goldstonc excitations in the 
stable phase where a quasiparticle gap arises. Let us show that the BS equation (|3.9[) indeed admits such solutions. 
To see this, according to 0], we set po = p = 0. Then, Eq. (|3.9[) has a solution of the form x(g, 0) = %5 (^0)75 for 
which we obtain the equation 



X5M) = ^ 



d 2 k 



x 5 (M) 



|q-k| A /k 2 + A 2 (fc)' 



(3.31) 



or, after integrating over the angle, 



X5(q,0) = A 



dk kxb(k, 0) 
v/fc 2 + A 2 (fc) 



K{q,k), 



(3.32) 



with the kernel 



K{q,k) = 



0(q-k) 



q 



9{k-q) ( q 



(3.33) 
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On the other hand, the equation for a gap function obtained in Ref. Q has the form 

A 




(3.34) 



One can see that the equation (|3.32[) has the solution X5(<?, 0) = CA(g) where the gap function A(q) satisfies the 
equation (|3 .34[) and C is a constant. Thus the wave function X5{Qi 0) describes a gapless Nambu-Goldstone excitation. 
Solving the BS equation at nonzero p 0l p one can obtain a dispersion law p ~ |p| for a Nambu-Goldstone excitation. 



In this paper we studied instabilities in graphene which arise at strong Coulomb coupling. For the supercritical 
Coulomb center problem, it was known before that the "fall into the center" instability arises if Za exceeds the critical 
value 1/2 leading to the appearance of quasistationary levels with complex energies. The energy of quasistationary 
states in the case of gapless quasiparticles has a characteristic essential-singularity type dependence on the coupling 
constant reflecting the scale invariance of the Coulomb potential. We showed that a quasiparticle gap stabilizes the 
system decreasing the imaginary part |Im£J| of quasistationary states, thus increasing their lifetime. 

Considering the many-body problem of strongly interacting gapless quasiparticles in graphene, we showed that the 
Bethe-Salpeter equation for an electron-hole bound state contains a tachyon in its spectrum in the supercritical regime 
a > a c and found the critical constant a c = 1.62 in the static random phase approximation. The tachyon states play 
the role of quasistationary states in the problem of the supercritical Coulomb center and lead to the rearrangement of 
the ground state and the formation of exciton condensate. Thus, there is a close relation between the two instabilities, 
in fact, the tachyon instability can be viewed as the field theory analogue of the "fall into the center" phenomenon 
and the critical coupling a c is an analogue of the critical coupling Z c a in the problem of the Coulomb center. The 
physics of two instabilities is related to strong Coulomb interaction. 

The calculated critical value a c — 1.62 should be compared with the value a c = 1.08 found in Monte Carlo 
simulations Q for the rearrangement of the ground state of graphene and appearance of a gap. The obtained value 
of a c is rather large that indicates that the ladder approximation is not quantitatively good enough for the problem 
of excitonic instability and gap generation in freely standing graphene. Certainly, both higher order corrections and 
improving the instantaneous approximation can vary the value of critical coupling. It is essential however that a ground 
state rearrangement at strong coupling is connected with the "fall into the supercritical Coulomb center" phenomenon. 
Therefore, such an rearrangement in graphene with large Coulomb interaction seems to be very plausible for strong 
enough coupling even if one goes beyond the ladder approximation. Finally, the phy sical picture of instabilities in 
graphene is quite similar to that elaborated earlier in strongly coupled QED [H, liH Il2| (see, also, [nj EH). In QED, 
the ladder approximation is not reliable quantitatively also because the critical coupling constant for chiral symmetry 
breaking is of order one. However, the main results of the ladder approximation survive when all diagrams with 
photons exchanges are included (the so-called quenched approximation without fermion loops) [3l|. Further, the 
existence of the critical point is exactly proved in the lattice version of QED [33|. We note also that in the presence 
of an external magnetic field the value of the critical coupling reduces to zero (magnetic catalysis phenomenon [35[) 
so that the gap generation takes place already in the weak coupling regime. 



We are grateful to V.A. Miransky and S. Sharapov for useful discussions. The work of E.V.G and V.P.G. was sup- 
ported partially by Ukrainian State Foundation for Fundamental Research under Grant No. F28. 2/083, by the grant 
10/07- N "Nanostructure systems, nanomaterials, nanotechnologies" , and by the Program of Fundamental Research 
of the Physics and Astronomy Division of the National Academy of Sciences of Ukraine. 

APPENDIX A: DISCRETE SPECTRUM FOR A REGULARIZED COULOMB POTENTIAL 

The discrete spectrum of Eq. ([2.4p exists for |e| < m. In this case it is convenient to define 
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p = 2ur, u 



V m 2 — e 2 , 



V m + e 



(<?-/), b 



ym — e 



(9 + /) 



(Al) 



a 



2 
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and rewrite Eqs. (|2.4p as follows: 



P9'+9\^-\- Za-) +f(j + Za-) = 0, 



2 2 uj V u 

[ P - + 1 --Za e -\ +g ( ] -Za™)=Q. (A2) 



2 2 u) * V u 

Substituting / from the first equation into the second one, we obtain the equation for the g component 

d 2 g (I h + Za t j-J 2 + Z 2 a 2 



dp 2 ' V 4 p p 2 

which is the well-known Whittaker equation [30,] . Its general solution is 



9 = 0, (A3) 



9 = C 1 W^[P]+C2M^[PI ^=1 + ^, (A4) 



where v = j 2 — Z 2 a 2 . Taking into account the asymptotic of the Whittaker functions Wk,u(z), Mk, v {z) at infinity. 

W„, v {p) ~ e~ ur (2ur)' i , r -> oo, (A5) 
M MtV (p) ~ -I(l±^e up (2«r)-", r ^ oo, (A6) 



we find that the regularity condition at infinity requires C2 = 0. Then the first equation in (|A2j) gives the following 
solution for the / component in the region II (r > R): 

hi = Ci (j - Za^) W_, +Za ^M (A7) 
Solutions in the region I (r < i?) can be easily obtained from Eqs 



°' = Al SgllC?) V itTJ/R-Z + f ) - ^ I ' (A9) 

where is a constant and we took into account the infrared boundary condition which selects only regular solution 
for bi and aj. Energy levels are determined through the continuity condition of the wave function at r = R, 

b -L\ r=R= b JL\ r ^ (AlO) 
ai an 



that gives the equation 

Wi + z*. (p) k + i m + e l e + Za/R-mJ\ j+1 /2\(p) „ m2 ^ 

t 5 — ^ — = , k = sen 7 1/ ; ■ ' „ , p = \/(Za + eR) z — m *R 

(j ~ Zf*) W_, + ^ v {p) r=R fc-l' 8 w u V e + Za/R + mJ^wip)^ Vl ; 

(All) 

We analyze this equation in the limit R — ► where we can use the asymptotical behavior of the Whittaker function 
at p — > 0, 

^00 - nl F(2 ^ + -**=^,H". (A12) 

In the limit R — > Eq. (|All[ ) reduces to the following one, 

■ Zalm+e) , ( ■ Za(m—e)\ 

T(-2v)T (l + v-Za^ ), - J + ^ V-^ + fc ob-^ J . , „, N 

" y (2ui?) 2,y = ) f + 0(i?), (A13) 



T{2v) T(l-v-Zaiy j v + k() 



Za(m — e) 
V s 
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FIG. 4: The lowest energy levels as functions of Za. Red lines correspond to the pure Coulomb potential (they exist for 
Za < 1/2); black solid lines are numerical solutions for j — 1/2, mR — 0.01; black dashed line are numerical solutions for 
j = -1/2, mR = 0.01. 



where 



Using the relationships 



m + e J\j + i/2\(Za) m + e 
sgn(j) J J lr7 <r(Za,j). (A14) 



u J\j-i/2\(Za) 



• . Za(m—e) rr • Za(m—e) r? 

j+v y - — '- Za u j - v y - — '- Za 



J u J 1 U 

can be rewritten in more convenient form 



(A15) 



r(2i/) r (1 - i/ - Za^) 1 ' j + u- Za ^- e ) j - v - ZaajZa, j) ' 1 ' 

In the limit R — ► the energy levels are determined by the poles of the gamma function r (l + v — Za^) and by a 
zero of the right hand side of Eq. (|A16[ ). this leads to the familiar result (analogue of the Balmer formula in QED) 
[lH (rederived also in (HI), 



1 



n = 12 1 3,...,i>0, (A1?) 
n=l,2,3,..., j<Q. y ' 



(u + n) 2 

The bound states for n > 1 are doubly degenerate, e nj = e n ,-j- The lowest energy level is given by 



eo,j=i/2 = my/1- (2Za) 2 . (A18) 

If Za exceeds 1/2, then the energy (|A18[) becomes imaginary, i.e., the fall into the center phenomenon (l8l. l23l l24j 
occurs. According to [ill. [l7|. nonzero R resolves this problem. For Za > 1/2, v is imaginary for certain j and for 
such j we denote v = i/3, j3 = \J Z 2 a 2 — j 2 . For finite R discrete levels exist for Za > 1/2. Their energy decreases 
with increasing of Za until they reach the lower continuum. The behavior of lowest energy levels with j = 1/2 as 
functions of the coupling Za is shown in Fig [4] 

The critical charge Z c that corresponds to diving into the continuum is obtained from Eq. (|Al6j) setting e = — m 
there and using the corollary of the Stirling formula: pfez|^ — * e 2tyloex , x — * +oo. We come at the equation 

- 2l /3 \og(2ZamR) = ifi - j + Za<r(Za,j) r(l - 2i(3) 

-i0-j + Zaur(Za,j)r(l + 2i0)' 1 ' 
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or, 

- (3 \og(2ZamR) = arg (Zaa(Za, j) - j + ij3) + arg T(l - 2i/3) + tto, (A20) 

where n is integer. It is not difficult to check that for j = 1/2 and n = 1 the critical coupling Z c a approaches the 
value 1/2 for mR — > 0. The dependence of the critical coupling Z c a on mR for j = 1/2 is shown in FigfT] 

The bound and quasistationary states in gapped graphene in the case of the supercritical Coulomb impurity were 
also numerically calculated in the tight-binding lattice model which has a natural lattice scale cutoff that provides an 
important control of the validity of the Dirac equation approach [U, [34[ . 



APPENDIX B: FOURTH ORDER DIFFERENTIAL EQUATION 



The system of equations (|3.17[) with boundary conditions (|3.18[) is reduced to the following fourth order differential 
equation for the function Xs{q)'- 



jy , 2(5 g 2 -3p 2 ) ,„ , (44 + 5A) 9 2 -8p 2 „ Ap 2 + (8 + 7X)q 2 , , 3A 2 X5 



q{q 2 -pl)' Xb ' 2q 2 {q 2 -p 2 ) 
with the corresponding boundary conditions 

d 



X 5 



q 3 (q 2 -pi) 



-x 5 



2q 2 (q 2 -p 2 ) 



0, 



<*Yb 

9=0 

(<7X5 (<?))' 



0, q" 



dq 



(q -pt) X5 + -X5 + A 



\5 



q 2 -Po 



9=0 



9=A 



0, 



dq 



q\q 2 ~pl) (x's' + ^ + A 



X5 



q 2 -Po 



9=A 



= 0. 



In terms of the variable z = q 2 jp\ these equations are rewritten as 



^ di „ 2/. on d 3 5 „„, „, d 2 19A + 60 d 3A 2 

z\ Z -l)- 4+ 2z 2 (4z - 3)^-3 + -z((X + 22), - 10)-^ + -^-z- + - 



X5 = 0, 



and boundary conditions 

,3/2 d X5 
dz 



= 0, 

z=o dz 



. d 2 x 5 . a dx 5 . Ax 5 



2z ~r~ + xi 

dz 



= 0, zV2" 
z=A 2 dz 



): 




= 0, 






2=0 


dX5 


Axs 


dz 


z-1 



z=A 2 



= 0. 



Eq. (|B4[) is the Pochhammcr-type equation [36|, its canonical form is 

/ 3 4 \ 



vfe=0 



fc=l 



where the parameters 



6 = 1, h = 3/2, b 2 = 3/2, 63 = 0, 



ai, 2 = -(l±7), 7 = Vl — 4A, a 3 ,4=-(l±7) 5 7 



2A 



(Bl) 

(B2) 
(B3) 

(B4) 

(B5) 
(B6) 

(B7) 

(B8) 
(B9) 



describe the behavior of X5( z ) & t the points z = and z = 00, respectively. The general solution of Eq. (|B4[) at the 
point z = can be written in terms of four linearly independent solutions, 



X.5 



Ci 



1-^3 Ol 



1 1 1 



> «2 — - , 03 — - , <X4 



2 2'2' 



, -, 1; z + C 2 Z4F3 ai + 1, a 2 + 1, a 3 + 1, a 4 + 1 



5 5 ^ 
; 2'2' ;Z , 



+ C 3 G 2 i ( z 



1 — 01, 1 — a 2 , 1 — a 3l 1 — a 4 
-1/2, -1/2, 0, 1 



C/^34 
4 ( - r 44 



1 — 01, 1 — 02) 1 — 03, 1 — a 4 
-1/2 0, 1, -1/2 / ' 



(B10) 
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where q +iF q ((a) q +i; (b) q ;z) is higher hypergeometric function and 



pq 



ai, . . . a n , a„+i, . . . a p 
b±, . . . b m , 6 m +i, . . . b q 



(Bll) 



is the Meijer G-function [32J. Leading asymptotic of the each term at z — > is 

4 



^=(1 + O(z)) + zC 2 (l + O(z)) + C 3 4=^- IT r ( a * ~ V2) (log ^ + £> + O(z)) 

+ a 4 ^-^nr(oi-l/2) + 0(l)Y D = 47-2 + 41og2 + ^V(ai-l/2), 
and 7 is the Euler constant. Hence from the boundary conditions (|B6|) we find that C3 = and 

"2 



X5 



d = C^Y[T(a^l/2). 



(B12) 



(B13) 



Asymptotical behavior of the function 4.F3 can be found from Eq.7.2.3.77 in the book [32j, thus we obtain 



n m) 



1-F3 (ai, a 2 , a 3 , a 4 ; bi,b 2 ,b 3 ; z) 



EH 



n r(oi) fc=i 



i=l 



r(a fe )nr«-a fe ) 

, t=i 

3 



(B14) 



if no two a*;, fc = 1, . . . 4, differ by an integer, the prime in the product Yi r( a '; — a,k) means that the term with i = k 

i=l 

is absent. Thus we obtain 



5 5 

Z4F3 ( 01 + l,a 2 + 1,03 + 1,04 + 1; -, 2' 2 ' 2 



T 2 (5/2) 
If r(a, + 1) 



(-*)" 



r(o 4 + l)r(oi - a4)r(a 2 - a 4 )r(a 3 - ai) 

r 2 (|-a 4 )r(i-a 4 ) 



(3 cyclic permutations 1^2^3^4^1) 



(B15) 



Similarly, for the Meijer G-function we use Eq.8.2.1.4 in the book [32j to find the asymptotic at large z: 



a±, . . . a n ; a n +i, . . . a p 
b±, . . . b m , & TO +i, . . . b q 



n%-<)nr(i+ii-a ( ) 
fe=i ft nak-bj) ft r(i+ aj - ak ) 



(Bi6) 



j—rn+l 



j=n+l 



if no two afc, fc = 1, . . . n, differ by an integer. 
In our case we obtain 



G 



.34 



1 — a>i, 1 — a 2 , 1 — as, 1 — a 4 
-1/2 0, 1, -1/2 



( z)~ ai r ^ ai - Q4 ^ r ^ 2 Z Q4 ) r ( Q 3 ~ a 4)r(a 4 )r(l + a 4 )r(a 4 - 1/2) 



r(3/2-a 4 ) 
(3 cyclic permutations 1^2^3^4-^1). 



Hence the function X5( z ) behaves at z — > 00 as 
4 

X5 (z) = A ^ (1 + 0(1/ z)) , Ai = (-l)- f - 



— C 4 7r 2 cot(7rai) — C 2 



V 



r 2 (5/2) 
ft r(a, + 1) 



Ft. 



(B17) 



(B18) 



where 



F, 



T(aj - a l )T(a k - aQT(oj - Qi)r(l + q^) 

r 2 (3/2-a 4 )r(i- ai ) 



(B19) 
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The UV boundary conditions lead to the following equations: 



- 7K 7/4 + A 2 (l + 7)^ 7/4 - A 3 (l + 3 7 K 1/2 ~ 37/4 - A 4 (l - S^z- 1 / 2 ^/ 4 = 0, (B20) 
At(l - 7 K 7/4 + A 2 (l + 7 )z 7/4 + A 3 ^(3 + 7 )z 1 /2-37/4 + ^3 (3 _ ^1/2+37/4 = 0j (B21) 

where z = A 2 /pq. This system of equations does not have solutions for A < 1/4. Near the critical value, A > 1/4, we 
can neglect the terms with A3, then we get A4 = and 



This gives 



and the equation 



Ai(l — ilo)z~ iuj I^ + A 2 (1 + iu)z tu ' A = 0, u = -if. (B22) 
r 2 (5/2) 

C 2 - V ' ' = -C 4 7r 2 cot(a 4 ), (B23) 

n r(o« + 1) 

i=l 



A A 1 — zcj i^i cot(7rai) — cot(-7ra4) 

^7 = 7 ( ) \- (B24) 

/ 1 + iu> F2 cot(7ra2j — cot(7ra4) 

From this equation we get when to <C 1, 

pi = -A 2 cxp[- — + a], aw 7.14, n=l,2,.... (B25) 

LU 

Comparing this result with Eq. (|3.24p we see that the approximation that decouples the system Q3.17P works nicely 
near the critical coupling A c . 
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